TSL/ISV-96-0152 



Scattering Wave Functions at Bound State Poles 



Goran FaldtQ 
Division of Nuclear Physics, Box 533, 751 21 Uppsala, Sweden 



Colin Wilking 

University College London, London, WCIE 6BT, UK 



February 1, 2008 
Abstract 

The normalisation relation between the bound and scattering S-state 
wave functions, extrapolated to the bound state pole, is derived from the 
Schrodinger equation. It is shown that, unlike previous work, the result 
does not depend on the details of the potential through the corresponding 
Jost function but is given uniquely in terms of the binding energy. The 
generalisations to higher partial waves and one-dimensional scattering are 
given. 
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1 Introduction 



It is demonstrated in standard quantum mechanics texts ||T[| that when a scatter- 
ing wave function is extrapolated to negative energy to the position of a bound 
state then the result is proportional to the bound state wave function. This is 
achieved through the divergence at this energy of the coefficient of the exponen- 
tially decreasing function at large distances. 

Much less interest has been placed on the relative normalisations of the scat- 
tering and bound state wave functions. Goldberger and Watson and subse- 
quently Joachain |^ have given a formal expression for this in terms of the Jost 
function and its derivative, which suggests that this quantity will in fact depend 
upon the form of the potential. 

We have recently shown [Q, using results from formal scattering theory, that 
the continuation of an S-state scattering wave function with real boundary con- 
ditions to a bound state pole actually only depends upon the bound state wave 
function and binding energy. This result is counter-intuitive in that a scatter- 
ing wave is normalised by its asymptotic behaviour at large distances, whereas a 
bound state is normalised by an integral condition. It is the purpose of the present 
work to give a self-contained demonstration of this starting from the Schrodinger 
equation and our principal result, eq. (|25|) , is to be found in §2. The generalisation 
of this to higher partial waves is outlined in §3, where the importance threshold 
kinematic factors in the extrapolation to the bound state is stressed. 

The relationship to the Goldberger and Watson form is investigated in §4, 
where it is shown that, carried to its logical conclusion, their form would give 
an answer identical to ours at the bound state position. Their extrapolation is, 
however, less smooth as a function of energy and, as a consequence, less useful in 
practice. 

The equivalent theorem for one-dimensional scattering is derived in §5 and 
some discussion of the usefulness of the results given in the conclusions of §6. 
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2 Derivation of the Theorem 



We start by summarising the relevant results of potential scattering as given, for 
example, in refs.[|, |]. 

Consider a particle of mass m moving with energy E = k!^ /2m and angular 
momentum £ in a real spherically symmetric potential V{r). In units where h = 1, 
the radial Schrodinger equation becomes 

^ Mk, r) + - Mk, r) = U{t) Mk, r) , (1) 

with U{r) = 2mV{r). 

In order that the results of standard scattering theory apply, we assume that 
the potential satisfies the integrability condition 

roc 

/ dr r I V{r) \ < oo . (2) 

^0 

This requires that the potential be less singular than at short distances but 
decrease more rapidly than at large ones. For such a regular potential the 
bound state eigenvalues of eq.(|Tp are non-degenerate and finite in number 

The generalisation to higher partial waves will be given in §3, and here we 
shall only consider S-wave solutions for which the index £ = on the radial wave 
function will be suppressed. 

The Jost solution f{k,r) is determined by the asymptotic behaviour 

f{k,r) ™ e""'. (3) 

With this boundary condition the radial Schrodinger equation (|I|) can be trans- 
formed into an integral equation and it is then straightforward to show by iter- 
ation that f{k,r) is analytic in the half-plane Im{fc} > and is continuous and 
bounded in Im{A;} > for all values of r > 0. 

The value of this solution at r = is the Jost function 

F(fc)^/(fc,0), (4) 
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and, unless this vanishes, the full wave function resulting from the Jost solution 
will be singular at the origin. 

In the scattering region, k is real and positive, whereas for a bound state 
k = ia, with a > in order that the bound state be normalisable. For any fixed 
value of r the scattering wave function f{k,r) can be analytically continued in 
k to give the bound state wave function f{ia,r) normalised such that it behaves 
like e~°'^ at large distances. 

A second independent solution of the Schrodinger equation is f{—k,r), and 
it is easy to see from the integral form of the equation that for Im{/c} > 

/(-r,r) = /(fc,r)*, (5) 

from which it follows that for the Jost function 

F{-k*) = F{ky. (6) 

For our purposes of continuing to the bound state pole, it is crucial to work 
with a real function. A physical radial wave function v{k, r) must vanish at the 
origin and the real linear combination with this property is 

v{k, r) = ^ [f{k, r) e''^'^ - fi-k, r) e'^^^'^)] , (7) 

with the phase shift 5{k) and corresponding S-matrix defined by 

S{k) = ^ = e-^(^') . (8) 

The function v{k, r) is real, due to the symmetry property of the Jost function 
shown in eq.(^, and behaves asymptotically hke 

v{k,r) — sin{kr + 6{k)) . (9) 
k 

The bound state positions are determined by demanding that the wave func- 
tion decrease at large r and be finite at the origin, and this is achieved if 
F{ia) = 0. The corresponding ^-matrix has a pole at this position and it is 
known that such poles are simple for a well-behaved potential satisfying eq.(0) 



For positive real values of k the norm of the S"- matrix is, by eq.(^, unity. We 
may parametrise it in the vicinity of an isolated bound state pole as 

S(k) ^ e^^^W = E^iMl! , (10) 
a + ik 

where N"^ is the residue at the pole. In the unitarised scattering length approx- 
imation N"^ = 2a and [(^(A;)]^ = {a — ik)/2a, but in general G{k) is an analytic 
function of k in the neighbourhood of A; = ia, with the condition that G{ia) = 1. 

Comparing the radial S-wave Schrodinger equations for the bound and scat- 
tering state wave functions, we have 

u"{r) — a'^u{r) = U{r)u{r) 
v"{k,r) + Pv{k,r) = U{r)v{k,r), (11) 

where U{r) is a real function. 

The boundary conditions at the origin are u{0) = v{k, 0) = 0, whereas the 
large-r behaviour of the bound state wave function is taken to be 

uir) e-"^ (12) 

while that of v{k,r) is given by eq.(^. The reality of the potential, combined 
with that of the boundary conditions, ensures that for real k the functions u{r) 
and v{k,r) remain real for all values of r. 

The wave function v{k,r) can be analytically continued in k to the bound 
state at k = ia, though care must be taken due to the singularity structure of 
the factor e^'^'-'^^ This is more complicated than its square, the S*- matrix, since it 
has a branch cut starting at the position of the bound state. 

It follows from eq. (p!OD that at the bound state pole 

J2a(a^ + k^)v{k,r)\ = - \Va + ik e'^^^^] , f(ia,r) = -Nu(r) . (13) 



k=ia 



It remains to be shown that the constant is determined uniquely by the nor- 
malisation of the bound state wave function. 
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Manipulation of the radial Schrodinger equations in eq.(]Tl|) leads to 

— [u'{r)v{k, r) — u{r)v'{k, r)] = (a^ + k'^)u(r)v(k, r) , (14) 

where prime indicates derivative with respect to r. Since the functions vanish at 
the origin, this can be integrated to give 

u'{r)v{k,r) — u{r)v'(k,r) = (a"^ + k'^) / dr'u{r')v{k,r') . (15) 

Jo 

Both sides of this equation vanish when r oo, reflecting the orthogonality 
of the scattering and bound state wave functions. To avoid a = statement, 
introduce the combination 

w{k, r) = 2ikV a + ikv{k, r) , (16) 

which, by eq.(0), has the limit at the pole 

w{ia,r) = N u{r) . (17) 

Differentiating the resulting 

u'{r)w{k,r) — u{r)w'{k,r) = (a"^ + k'^) / dr'u{r')w{k,r') (18) 

Jo 

with respect to k leads to 
u{r)w{k,r)—u{r)w'{k,r)= / dr'u(r') {o^ + k^)w{k,r') + 2kw{k,r') , (19) 



where derivatives with respect to k have been indicated by dots. 

After taking first the limit k —>■ ia and then r —>■ oo, the first term in the 
integrand vanishes because of the explicit (a^ + k'^) factor, and so the right hand 
side of eq . (|I9|) becomes 

POO 

- 2iNa / dru'^{r) . (20) 



On the other hand, from the ansatz of eq. ([T0|) , in the vicinity of the pole 

w{k,r) = NG{k)f{k,r)~^^±^f{-k,r), 

w{k,r) = NG{k)f{k,r) + NG{k)f{k,r)-j;^^^f{-k,r) 

fi-k,r) + ^f^T7TX^G{k) f{-k,r) . (21) 



NG{k) ' ' NG{kY 
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Thus, since G{ia) = 1 at the pole, 

w{ia,r) = Nf{ia,r) ™ A^e""", (22) 
w{ia, r) = N G{ia) f{ia, r) + N f{ia, ^) ~ f{~i<^^ ^) 

™ _i_e+-. (23) 

Since the bound state wave function decreases exponentially, the only surviv- 
ing contributions to the left hand side of eg . ([191) at large r originate from the 
term proportional to f{—ia,r) in eg. (p3|) , yielding a value of —2ia/N. 

Eguating the two sides of eg.(|TP|), leads to the condition that 

1 = N^ dru\r), (24) 





so that is indeed the normalisation constant of the bound state wave function. 

Rewriting the result in terms of a normalised radial bound state wave function 
Ua{r) = Nu{r), we see that the scattering wave function v{k, r) is related to Uair) 
at the pole k = ia through 

lim lj2a{a^ + k^)v{k, r)\ = -M«(r) , (25) 

and this is the principal result of this work. 

It should be noted that, in contrast to the result given in ref. Q], the constant 
of proportionality in this relation is independent of the form of the potential 
and of the Jost function. Note also that, since a particular bound state is not 
specified, eg.(p5|) is valid at all bound state poles. 
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3 Higher partial waves 



For higher ^-values the discussion given in §2 can be repeated with only minor 
changes. The Jost solutions of eq. (|^) are determined by the asymptotic behaviour 

i 

fi{±k,r) Qini/2^±ikr ^ (26) 

and the Jost functions defined as 

The two Jost solutions are related for real k by 

f,(-k,r) = {-lY[Mk,r)r , (28) 

and similarly for the Jost functions. 

The physical radial wave function ve{k, r) must vanish at the origin like r^"*"^, 
and the real linear combination with this property is 

ve{k, r) = ^ [{-lYfiik, r)e^^^('=) - ^-k, r)e-^^^(^)] , (29) 



which for large values of r tends to 



Vi{k,r) ™ l-sm{kr-li7i + 6e{k)) . (30) 



In the vicinity of the bound state pole we write the S-matrix 

FfXk) a + ik 



SKl0^e-M'. = fMzM=(zim££W]!, (31) 



with Ge{ia) = 1 and A'^ real. 

The residue has an ^-dependent phase factor and we define the square root 
through 

ya + ik 



The extra phase factor matches that in the definitions of eqs. f^Bf ) and (p^) 
such that, with Wi{k,r) defined as in eq. ([T6|) , the results of eqs.([17|) and (pT^) 
remain unchanged with Ni the asymptotic normalisation constant for angular 
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momentum In particular, the bound state wave function with the asymptotic 
normahsation of eq. ([T^) is defined by 

ui,a{r) = {-iYfi{ia,r) . (33) 

In sum, our principal result of eq. remains true for arbitrary partial waves. 
However, in order to provide a smooth extrapolation to the bound state pole it is 
preferable to take out the explicit threshold behaviour by dividing by a factor 
to give 

^lim |v'2«(«2 + p) (^|y^^(fc^^)| = _(_,)^^^^^(^) . (34) 

Though of little importance in practical applications, this introduces an extra 
overall phase factor, whose origins may be traced to the definition of the partial 
wave solution of eq. (pQ]) . 
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4 Relation to the Goldberger and Watson Form 

Using formal scattering theory arguments, Goldberger and Watson have de- 
rived an expression relating the normalisation of scattering and bound state wave 
functions |@| for the regular S-wave solution ip{k,r), which is fixed by two real 
boundary conditions at the origin. The result, which depends upon the values of 
the Jost function and its derivative at the S-matrix pole, cannot be simplified, 
which means that the relation between the regular solution ip{k, r) and the bound 
state wave function Ua{r) depends intrinsically upon the form of the potential. 

The Goldberger- Watson form can however easily be recast to give a relation 
for the physical scattering wave function v{k,r) {cf. eq.(11.89) of ref.0). 



implementation of eg . (|36|) . 

However for real values of k our expression for the extrapolation function in 
eq. (l25|) remains explicitly real, whereas that of eq.(^) becomes complex. As a 
consequence the former provides a smoother extrapolation for the relationship 
between two real quantities. 

To see the difference in practice, consider a one-term separable potential with 
shape parameter f3, for which the ratio of the scattering and bound state wave 




(35) 



In the vicinity of an isolated zero the Jost function may be written as 



F{k) =iC{k) (k-ia), 



(36) 



where C{k) is regular and C{ia) ^ 0. 
It then follows immediately that 
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functions at the origin is given by 



i^W-l^l (38) 



where 

A;2 + ^2 

^ = 4/3(a + /3) ■ ^^^^ 
The correction terms to eq.(^) are therefore small at r = providing /c^ and 
are small compared to 

The extrapolation function for the Goldberger and Watson form may be easily 
derived from the Jost function 

^ {k-iP)[k + i{2^ + a)]{k-ia) 
^ ' {k + il3W + {a + l3Y + l3^] ■ ^ ' 

It is easily seen that this leads to large errors as soon as k'^ is comparable to 
and this is due to the neglect of the influence of the zero of the S-matrix at 
k = —ia. 
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5 One-dimensional scattering 



One-dimensional scattering problems show certain features which are not present 
in the more usual three-dimensional case 0, but our proof goes through here 
largely unaltered. 

Assuming the potential to be an even function of the variable x, for which 
V{—x) = V{x), it is convenient to work with even and odd solutions of the 
Schrodinger equation v^{k, x) and v^{k, x) which at the origin satisfy the bound- 
ary conditions 

<(fc,0) = 0, (41) 
v_{k,0) = 0. (42) 

For large positive values of x the functions behave like 

v+{k,x) ^-^^ cos{kx + 6-^-{k)) , (43) 
v_{k,x) ™ -sm{kx + 6^{k)) . (44) 

We define the normalisations of the corresponding symmetric and antisym- 
metric bound state wave functions through 

(45) 

With the exception of an overall —k factor in eq.(0), the boundary condi- 
tions on the antisymmetric function are identical to those of the S-wave three- 
dimensional problem of eq.(|^) combined with the vanishing at the origin. Taking 
into account that the normalisation integral extends down to — oo in eq. (^5]) , one 
can deduce the one-dimensional extrapolation directly from the corresponding 
three-dimensional case in eq. (|25|) 

lim y'a(a2 + F) v{k, x)] = m"(x) . (46) 

The presence of the 1/k factor in eq.(^6D does not upset the smoothness of 
the extrapolation since this is compensated by the vanishing of v^{k, x) at k = 0, 
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r+oo 




2 r+oo 




2 


/ dx 


ulix) 


= dx 


uHx) 


= 1 


/ — oo 




J —oo 







which is a clear consequence of the antisymmetry of the function with respect to 
kx. 

For the symmetric case we must modify the arguments given in §2. In the 
vicinity of a simple pole of the S-matrix, we may write 

= (47) 

(y. ~\~ trv 

where G{k) is analytic in the neighbourhood of the pole and is normalised there 
to G{ia) = 1. 

Now it is known that as A; ^ 0, S^{k) 7t/2 unless there is a zero-energy 
bound state . This is one of the special features of one-dimensional scattering, 
which causes the one- dimensional transmission coefficient to vanish at zero energy. 
When this condition is imposed upon eg. (^7D , it follows that 

G(0) = — ,/^- (48) 

Due to the symmetry of the problem, we need only consider the region a; > 0, 
and there the Jost solution is defined by its asymptotic behaviour 

f{k,x) ™ e'^". (49) 

This allows the even function to be written as 

v{k, x) = ^ [/(A;, x) e'^+(^) + /(-A;, x) e-'^+(^)] . (50) 

After defining a non-normalised bound state wave function m_|_(x) with a large- 
X limit of e~°^, as in eg. (P^) , all the subseguent manipulations follow as in the 
three-dimensional proof. This is because in the symmetric case the derivatives of 
the wave functions vanish at the origin, and so it is still permissable to discard 
the left hand side of eg. ([T5|) at x = to leave, as before, 

rx 

u'_^_{x)v^{k,x) — u^{x)v'_^{k,x) = {a'^ + k'^) / dx'u^{x')v^{k,x') . (51) 







It is then straightforward to show that the one- dimensional result of eg. 
is valid for both symmetric and antisymmetric wave functions. 
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6 Conclusions 



We have given a simple proof of the relation between the normalisations of the 
bound and scattering wave functions at the bound state pole which is on a par 
with the standard proof of the effective range expansion [|l|. Apart from being 
formally correct, the result is useful in practice because the extrapolation function 
is smooth in k"^. As we saw in the case of the separable potential in §3, corrections 
to it only become important when 1/ A; is of the order of the range of the potential. 
It is therefore possible to use the theorem to form approximations for scattering 
wave functions in terms of that of the bound state. For the S-wave Paris neutron- 
proton triplet potential such an approach is reliable out to radii of 1.7 fm for 
centre-of-mass energies up to at least 20 MeV [§. 

Already in 1952, Watson used such arguments to estimate the final state 
interaction effects in pp n^i^pn) at low excitation energies in terms of the 
cross section for pp — > vr+rf. The relative scattering/bound state normalisation 
was there established as an approximation using effective range theory, though 
in practice it differs little from our result at the pole. Providing the momentum 
transfer is large, as it is for pion production, then the corresponding transition 
operator is sensitive to the short range part of the neutron-proton wave function 
where our theorem provides a valid approximation. In addition to providing a 
useful description of pp -K^ipn) data near threshold and at higher energies, 
it allows one to understand quantitatively the final state interaction regions in 
pn — > ri{pn), pd — > TT^{pd), and dp —>■ p{pn) at high momentum transfers @, 10]. 

The extension to higher partial waves described in §3 is likely to be of less 
significance since it is harder to investigate the effects of P-wave final state in- 
teractions unless the S-wave is suppressed. Nevertheless the modifications to our 



principal result of eq.(p5D, engendered by threshold kinematic factors for £ 7^ 0, 
are of importance here. 
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